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ABSTRACT 

We re-examine the outer gap size by taking the geometry of dipole magnetic field into account. 
Furthermore, we also consider that instead of taking the gap size at half of the light cylinder radius 
to represent the entire outer gap it is more appropriate to average the entire outer gap size over the 
distance. When these two factors are considered, the derived outer gap size f(P,B, < r > (a)) 
is not only the function of period (P) and magnetic field (B) of the neutron star, but also 
the function of the average radial distance to the neutron star < r >; which depends on the 
magnetic inclination angle (a). We use this new outer gap model to study 7-ray luminosity of 
pulsars, which is given by L 1 — f 3 (P, B,< r > (a))L s d and L s d is the pulsar spin-down power, 
as well as the death lines of 7-ray emission of the pulsars. Our model can predict the 7-ray 
luminosity of individual pulsar if its P, B and a are known. Since different pulsars have different 
a, this explains why some 7-ray pulsars have very similar P and B but have very different 7-ray 
luminosities. In determining the death line of 7-ray pulsars, we have used a new criterion based 
on concrete physical reason, i.e. the fractional size of outer gap at the null charge surface for 
a given pulsar cannot be larger than unity. In estimate of the fractional size of the outer gap, 
two possible X-ray fields are considered: (i) X-rays are produced by the neutron star cooling and 
polar cap heating, and (ii)X-rays are produced by the bombardment of the relativistic particles 
from the outer gap on the stellar surface (the outer gap is called as a self-sustained outer gap). 
Since it is very difficult to measure a in general, we use a Monte Carlo method to simulate 
the properties of 7-ray pulsars in our galaxy. We find that this new outer gap model predicts 
many more weak 7-ray pulsars, which have typical age between 0.3-3 million years old, than the 
old model. For all simulated 7-ray pulsars with self-sustained outer gaps, 7-ray luminosity L 7 
satisfies L 1 cx L s sd ; where the value of 5 depends on the sensitivity of the 7-ray detector. For the 
EGRET, 6 is ~ 0.38 whereas S is ~ 0.46 for the GLAST. For 7-ray pulsars with L sd < Lf d lt , 
S is ~ 1. Lf d lt = 1.5 x 10 34 P 1 / 3 ergs _1 is determined by /(< r >^ tl) = 1. These results 
are roughly consistent with the observed luminosity of 7-ray pulsars. These predictions are very 
different from those predicted by previous outer gap model, which predicts a very flat relation 
between L 1 and L s d- 

Subject headings: pulsars: general-star: neutron- gamma-rays: stars 

1. Introduction els, charged particles are accelerated in charge- 

depleted zones near the pulsar's polar cap and 
High-energy emission models for rotation- 7 _ rays are produced through curvature-radiation 

powered pulsars are generally divided into polar induced 7 -B pair cascade (e.g. Harding 1981; 

gap and outer gap models. In polar gap mod- Daugherty & Harding 1996; Zhang & Harding 



1 



2000) or through Compton-induced pair cascades 
(Dermer & Sturner 1994). In the outer gap mod- 
els, it is generally accepted that a magnetosphere 
of charge density 



surrounds a rotating neutron star with magnetic 
field B and angular velocity SI (Goldreich & Julian 
1969). The magnetospheric plasma is corotating 
with the neutron star within the light cylinder, at 
which the corotating speed equals the velocity of 
light and the distance from the spin axis is 

Rl = c/fi. In the corotating magnetosphere, 
the electric field along the magnetic field, E\\ = 
E • B/_B, is nearly zero. However, the flows of the 
plasma along open field lines will results in some 
plasma void regions (where the charge density is 
different significantly from p ) in the vicinity of 
null charge surfaces where ft ■ B — (Holloway 
1973). In such charge deficient regions, which 
are called outer gaps, E\\ ^ is sustained, elec- 
trons/positrons can be accelerated to relativistic 
energies and the subsequent high-energy gamma- 
ray emission and photon-photon pair production 
can maintain the current flow in the magneto- 
sphere (Cheng, Rudcrman & Sutherland 1976; 
Cheng, Ho & Ruderman 1986a, 1986b, hereafters 
CHRI and CHR II; Romani, 1996; Zhang & Cheng 
1997; Hirotani 2001 ). 

Based on known 7-ray pulsars, the luminosity 
and conversion efficiency of 7-rays in various mod- 
els have been studied (for example, Harding 1981; 
Dermer & Sturner 1994; Rudak & Dyks 1998; 
Yadigaroglu & Romani 1995; Zhang & Cheng 
1998). Observations by Compton Gamma-ray Ob- 
servatory (CGRO) show that 7-ray luminosity of 
rotation-powered pulsars is proportional to square 
root of the spin-down power (Thompson et al. 
2001). Using the recent new polar gap models 
(e.g. Zhang & Harding 2000; Harding & Mus- 
limov 2001; Harding, Muslimov & Zhang, 2002), 
Harding et al. (2002) have studied the deadline 
of 7-ray pulsars based on the predicted luminosity 
of 7-ray pulsars L~ ( cx L s sd , where 8 <~ 0.5 when 
Lsd > L b / d eak and 5 ~ 1. when L sd < L b / d eak re- 
spectively and L b ; d eak = 5 x 10 33 p-V2 erg / s . 

For a rapid rotating pulsar, it is believed that 
its spin-down power, L s d, is converted into radia- 
tion energy. Because the outer gap occupies only a 



part of the open field line region, the gamma-ray 
luminosity produced in the outer gap is a frac- 
tion of the spin-down power. It has been shown 
that the gamma-ray luminosity in the outer gap 
is proportional to / 3 , i.e L 1 w f 3 L sc i (CHR II; 
Zhang & Cheng 1997). In previous works, the fac- 
tional size of outer gap only depends on the pe- 
riod and magnetic field for a gamma-ray pulsar. 
For example, the fractional size of the outer gap 
for Crab-like pulsars is / cx B~ 2 13/2 ° P 33 / 20 (CHR 
II). In the outer gap model described by Zhang 
& Cheng (1997), / cx B^/ 7 P 26 / 21 , the observed 
gamma-ray luminosities with energies greater than 
100 MeV from the known gamma-ray pulsars ex- 
cept for the Crab pulsar can be explained approxi- 
mately (Zhang & Cheng 1998). It should be noted 
that the model predicts that the 7-ray pulsars with 
the same values of (B/P) - 3 have same 7-ray lu- 
minosities. For example, the ratio of (B/P) 3 for 
PSR B1055-52 to that for Geminga is - 0.9, it 
means that the ratio of both gamma-ray luminosi- 
ties is ~ 0.9. However, the observed ratio of the 
luminosities with energies greater than 100 MeV 
of these two pulsars are ~ 8 (Kaspi et al. 2000). 
Briefly, 7 known 7-ray pulsars have rather different 
7-ray luminosity even their spin-down powers and 
ages are so similar (e.g. Geminga and PSR 1055- 
52). Their pulse shapes also differ so much. It 
is clear that there are other intrinsic parameters 
to control these observed properties (gamma-ray 
luminosity, pulse shape , spectrum etc). 

In this paper, we re-study the gamma-ray emis- 
sion from the outer gaps of the rotation-powered 
pulsars by using a new outer gap model. We follow 
the idea of self-sustained outer gap given by Zhang 
& Cheng (1997). However, we take the magne- 
tosphere geometry as well as the average proper- 
ties of the entire outer gap into consideration and 
show that the fractional size of the outer gap is 
a function of period, magnetic field and magnetic 
inclination angle. In fact, the effect of the incli- 
nation angle on the 7-ray emission have been con- 
sidered in other versions of the outer gap models. 
For example, Romani & Yadigaroglu (1995) and 
Yadigaroglu & Romani (1995) took the magnetic 
inclination angle into account in their outer gap 
models, Hirotani and his colleagues also included 
the magnetic inclination angle in their calculation 
of the outer gap model (e.g. Hirotani, 2001; Hi- 
rotani & Shibata 2001; Hirotani & Shibata 2002: 
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Hirotani, Harding, Shibata 2003). Different from 
the treatment of Romani & Yadigaroglu (1995) 
who considered it in a less analytic way, we give 
a explicit expression for the fractional size of the 
outer gap. In section 2, we describe the revised 
outer gap model. We estimate the 7-ray lumi- 
nosity for rotation-powered pulsars and compare 
them with the observed data in section 3. In sec- 
tion 4, we derive the death lines of the pulsars with 
outer gaps. Finally, we give briefly our conclusion 
and discussion. 

2. The Outer Gap Model 

2.1. Magnetospheric Geometry 

For an oblique magnetic dipole rotator with an 
angular velocity f2 and the magnetic moment vec- 
tor yit, let its spin axis be along the Oz axis, \x be 
in the plane xOz and a be the angle between fi 
and \i. In polar coordinates, 



n = Q(cos0f -sin 00) . 



(2) 



and \i = /u(cos(0— a)f, — sin(0 — a)0). Correspond- 
ing magnetic field is 

B(r) = ^g(2cos(0-a)r + sin(9 - a)9) , (3) 

where \i = B p R 3 /2, B p and R are the stellar ra- 
dius and surface magnetic field at pole (see, for 
example, Zhang & Harding 2000) . f and are the 
unit vectors of radial and polar angle directions 
respectively. For a pulsar with a period P and a 
period derivative P, the B p is estimated by 



B p w 6.4 x 10 19 (PP) 1/2 G . 



(4) 



It is believed that the outer gap is extended 
from its inner boundary to the light cylinder (CHR 
I). For the oblique magnetic dipole rotator, the 
polar angle, C , in which the last open field line is 
tangent to the light cylinder given by (Kapoor & 
Shukre 1998) 



tan0 c = — - 



4 tan a 
corresponding radius is 



(l + (l + 8tan 2 a/9) 1/2 ) , (5) 



Rl 
sin 9 r 



(6) 



it should be noted that 9 C = n/2 and r c = Rl for 
an aligned magnetic dipole. Along the last open 
field line, the relation 



sin 2 (0 — a) sin 2 (0 c — a) 



(7) 



is valid. The inner boundary of the outer gap is 
estimated by the null charge surface, which is de- 
fined by $7 • B = 0. In the two dimensional case, 
the null charge surface can be described by (r irl , 
9 in ). By definition, we have 



tan 6 i n = — (3 tan a + \J 9 tan 2 



(8) 



while ri n is estimated along the last open field line, 
which gives using Eq. (7) 



rin _ sm 2 (9 in - a) 
Rl sin 9 C sin (0 C — a) 



(9) 



Therefore, the outer gap extended from n n to r c 
along the last open field line for the oblique mag- 
netic dipole. In such a geometry, the Goldrcich- 
Julian current is roughly 



<t ~ Q 2 R 3 B P 

— 2ec — °° Sa 

where a(a) — sin0 c sin 2 (0 c — a). 



(10) 



2.2. X-ray Field in the Magnetosphere 

The observed spectra of X-ray emission from 
some rotation-powered pulsars indicate that there 
are at least two kinds of X-ray spectra. One con- 
sists of soft X-rays, which can be fitted by a black- 
body spectrum with a single temperature, com- 
bined with a hard X-ray spectrum with a power- 
law distribution. Another has only a thermal spec- 
trum. It is believed that the non-thermal com- 
ponents are most likely of magnetospheric origin, 
while the origin of thermal components is less clear 
because there are several mechanisms for produc- 
tion of thermal emission in the soft X-ray bands. 
Observationally, the bulk of the soft X-ray emis- 
sion between <~ 0.1- 1 kcV is well fitted by a dou- 
ble black-body at two different temperatures. The 
ratio of the area of the hotter component to the 
colder one is typically very small (~ a few xl0~ 5 
to a few xl0~ 5 ). The colder component has been 
explained as resulting from thermal cooling, while 
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the hotter one most likely comes from bombard- 
ment by high-energy particles (Greiveldinger et al. 
1996). It is expected that pulsars younger than 
~ 10 5 yr have significant neutron star cooling com- 
ponents. We describe briefly the main possible 
mechanisms of X-ray production as below. 

2.2.1. Thermal X-ray emission from neutron 
star cooling 

For thermal X-ray emission due to neutron star 
cooling, it is believed that its spectrum is ex- 
pressed with a modified black-body, we approxi- 
mate it as a black-body spectrum with a temper- 
ature T c . Because the neutron star cooling con- 
cerns many different mechanisms, the estimate of 
the temperature T c exists uncertainties for differ- 
ent models. For example, the temperature can be 
expressed as (Romani 1996) 



'c,6 



10 5 yr 



-0.05 



exp(-T/10V) K , (11) 



where T C:6 is the temperature of the neutron star 
cooling in units of 10 6 K and r is the neutron 
star's age. Above expression is valid for r less 
than several 10 6 yr. Zhang & Harding (2000) used 
a expression to approximate the temperature T c as 

_ f 10-°- 23 r 6 - 01 r<10 5 V n9] 

Jc,6- | 10 -0.55 T -0.5 T>10 5.2 yr; W 

where tq is the age in units of 10 6 yr. Based on the 
cooling model derived from Tsuruta (1998) (in this 
model, the effects of the magnetic field on thermal 
conductivity are included, but polar cap heating 
is not considered), Hibschman & Arons (2001) use 
a temperature model as follow 



TcAt) 



10 
10 



0.05„-0.1 
0.325--0.375 



r < 10 7 yr 
t > 10 7 yr 



(13) 



where t is the spin-down age of the neutron star. 
For the hotter component originating from the 
bombardment by high-energy particles, the high- 
energy particles are produced either in the polar 
cap region or in outer gap region, resulting in 
different properties of the hotter component. In 
other words, the origin of the relativistic parti- 
cles, which bombard the stellar surface to produce 
thermal hotter X-rays, depends on the models. 



2.2.2. X-ray emission from polar cap heating 

In the polar cap models, the return relativistic 
particles are produced in the polar gap (e.g. Rud- 
erman & Sutherland 1975; Arons 1981). Recently, 
Zhang & Harding (2000) considered full polar cap 
cascade scenario of the polar cap model and esti- 
mated thermal 

X-ray luminosity using a self-consistent polar 
cap heating in the Harding & Muslimov (1998) 
model (a relevant recent study see Hibschman 
& Arons 2001). Further, Harding & Muslimov 
(2001) studied the effect of pulsar polar cap heat- 
ing produced by 

positron returning from upper pair formation 
front, the polar cap heating produces a ther- 
mal X -Tcty emission with ci temperature, Tp C — 
(Lpc/aA) 1 ^, 



2.46 



(Pn i A ™ ( cos 2 a 
\ t 6 J \ a(a) 



forP^ < 0.5^12 
forPj, > 0.5^12 



(14) 

where L pc is the X-ray luminosity emitted from 
the polar caps, a is the Stefan constant, and A 
is the heated area of the returning positrons from 
the polar gaps. For the canonical polar cap area, 
A = A pc = irR 2 (nR/c). In the above expressions, 
P 01 = P/0.1 s , B 12 = B/10 12 G and t 6 = 10 6 r 
yr is the pulsar age, we have used the analytic es- 
timate of X-ray luminosity given by Harding and 
Muslimov (2001) and A pci so the above expres- 
sions are valid for normal pulsars with r < 10 7 
yr. 

2.2.3. X-ray emission from outer gap heating 

In the outer gap models, part of the relativistic 
particles from the outer gap will collide the stel- 
lar surface, producing the thermal X-rays. These 
relativistic inflowing particles from the outer gap 
radiate away much of their energy before reaching 
the polar cap (Zhang & Cheng 1997, Zhu et al. 
1997; Wang et al. 1998; Cheng & Zhang 1999). 
The residual energy of the charged particles strik- 
ing the polar cap is 



E e (R) 



2e 2 c r 

5 ' n T: 

m&Rh R 



-1/3 



(15) 



These particles collide with the polar cap at a rate 
of iV e = fNaj, where Nqj is the Goldreich- Julian 
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particle flux (Goldreich & Julian 1969) and es- 
timated by Eq. (10). Therefore the polar cap 
is heated and radiates X-rays with a luminosity 
(L x w fE e (R)N GJ ) 

L x « 2.3 x 10 31 /-Bi 2 P" 5/3 

i? 6 ^ln— J a(a) cos a ergs s _1 (16) 

These X-rays have a typical temperature Th = 
{Lx/Aass) 1 ^ with 

T h « 5.0 x lO^- 1 / 6 ^ 4 ^-^)" 1712 

(sin 2 (a(a) cos a) 1/4 K , (17) 

where A <~ nfClR 3 /c is the area of the stellar sur- 
face bombarded by the return current and ctsb 
is Stefan's constant. Generally, a part of these 
X-rays will escape along the open magnetic field 
lines. Following Cheng & Zhang (1999), we use 
L x as the luminosity of the X-rays escaping along 
the open field lines and introduce a parameter 
£ = L x I Lx (the estimate of £ see Cheng & Zhang 
(1999)). Most of X-rays are reflected back to the 
stellar surface because of the cyclotron resonant 
X-ray reflecting mirror. This process transfers 
emitted polar cap X-ray energy to the entire sur- 
face of the neutron star. Finally, the X-rays from 
the entire surface of the neutron star are emit- 
ted with a temperature of T = (L x /2ttR 2 o-) 1 / 4: , 
where L s x = (1 — £)Lx- Corresponding charac- 
teristic temperature is T s — (L x /4:iTR 2 asB) 1 ^ 4: , 
which gives 

l/4 f l/4p-5/12 R l/4 



4.2 x 10 b (l - 7 / 1 P 

-1/12 
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(ln-^) (a(a)cosa) 1/4 i?g /4 K(18) 

2.2.4- Average energy of X-rays 

We consider two possible cases for the thermal 
X-rays from the stellar surface. In the first case, 
the thermal X-ray are produced by both the neu- 
tron star standard cooling mechanism and polar 
cap heating, the temperatures are given by Eq. 
(12) or Eq. (13) for the standard cooling mecha- 
nism and Eq. (14) for the polar cap heating. In 
the second case, the thermal X-ray come from the 
bombardment of the relativistic particles from the 
outer gap (e.g. Zhang & Cheng 1997), the cor- 
responding temperatures are given by Eqs. (18) 



and (17). Assuming these X-ray can be approx- 
imated as the black-body, their spectrum can be 
expressed as 



Fx(Ex) — C 



+ 



E\ 



e E x lkT x _ l 4lrR 2 e E x /kT 2 _ l 

(19) 

where R is the stellar radius, A is the area of polar 
cap being heated, and T\ and T 2 are the temper- 
atures being whole stellar surface and polar cap 
respectively. Using above expression, we can esti- 
mate the average X-ray energy as 



< E x > 



J Fx{E x )E x dEx _ 7T 4 



jFx(E x )dE x 30C(3) 



1 + {A/A-kR 2 )(T 2 /T{)- 
1 + (^/47ri2 2 )(T 2 /Ti) : 



;(fcT 1 )(20) 



where £(x) is the Zeta function and £(3) w 1.2. 
For the different mechanisms of thermal X-rays 
from the stellar surface, the dependence of < 
Ex > on the basic parameters of pulsars is dif- 
ferent. 

We consider two possible cases for estimating 
the average X-ray energy: (i)thermal X-rays are 
produced by the standard neutron star cooling and 
polar cap heating, we have 



< X > ~ 2Jklc 1 + {R/AR L ){T PC /T c f 



(21) 



and (ii)the thermal X-rays are produced by the 
bombardment of the relativistic particles from the 
outer gap, we have 



< E° 9 2.7kT s ^-4: 
x 1 - £ 



(22) 



2.3. The Fractional Size of an Outer Gap 

In two dimensional geometry, the fractional size 
of the outer gap is an important parameter for the 
7-ray production in the outer gap. According to 
Zhang & Cheng (1997), the parallel electric field 
in the outer gap can be approximated as 



E\\ - fB(r) 



s 

Rl 



(23) 



where / is the fractional size of the outer gap, 
B(r) = B p (l + 3cos 2 (6 -a^/iR^/r 3 is the mag- 
netic field strength at the radius, r, to the star, 
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Rl is the radius of the light cylinder and s is the 
curvature radius which is (Lesch et al. 1998) 



8(0,0,) 



R sin(0-a) (1 + 3 cos 2 (<9 - a))i 
3"sin 2 (6» s - a) I + cos 2 (6 -a) 



(24) 

where 9 S is the polar angle at the stellar surface. 
Equation (25) can be written by 



8(0,0.) = y^Rl\V(a,r) 



with 



W(a, r) = — - 



3 yfi@j(l -&(<*)£) 



(25) 



(26) 



where a(a) — sin 2 (9 c — a) sin# c . This electric field 
will accelerate the electrons/positrons to relativis- 
tic energy in the outer gap. Because these accel- 
erated particles will lose their energy through cur- 
vature radiation, their Lorentz factor is estimated 
by using eE^c = (2/3)e 2 c-f 4 / s 2 , which gives 



7(r) 



2.84 x 10 7 fiB? 2 P-iRl (^-j 



^/l + 3cos 2 (9 -a) S 



(27) 



where P is the pulsar period in units of second 
and R is the stellar radius in units of 10 6 cm. The 
characteristic energy of the 7-ray photons in the 
outer gap can be approximated as 



£ 7 « 143/^P"* 



vT + 3cos 2 ((9-a) 



r \ 2 / s V 



R'l MeV , (28) 



Inside the outer gap, the curvature photons 
interact with the thermal X-rays from the stel- 
lar surface to produce pairs through photon- 
photon pair production process, sustaining the 
outer gap. This pair production condition is 



2\2 



< E x > £ 7 (1 - cos(6> X7 )) = 2(m e cT) 



(29) 



where < Ex > is the average X-ray energy which 
is estimated below, 9xj is the angle between the 
emission directions of curvature photons and the 



thermal X-rays. We assume that the curvature 
photons are emitted along the negative direction 
of the magnetic field and the thermal X-rays along 
the radial direction, then we have 



cos 9x^ 



2cos(# — a) 



(3cos 2 (0-a) + l) 1 /2 ' 



(30) 



where 9 is the polar angle at radius r. Putting Eq. 
(7) into Eq.(30), we have 



cos 9 



x 7 



1 - (r/R L )a(a) 
(1 - (3/4)(r/R L )a(a)) 



1/2 



(31) 

where a(a) — sin (9 C — a) sin# c . 

We consider the fractional sizes of the outer gap 
corresponding to two possible average energies of 
X-rays. In the first case, X-rays are produced by 
the neutron star cooling and polar cap heating. 
The average X-ray energy is given by Eq. (21), 
and the fractional size of the outer gap is 



f(r,a) w 6.9B 12 *Pe < E x > .F G pc (r,a) 



(32) 



with 



G pc = W 1 ^ 



COS I 



X7 



2/3 



v / l + 3cos 2 (9 - a) 



(33) 



where < Ex >o.i=< Ex > /0.1 keV. Because the 
temperature of the polar cap heating is greater 
than that of the neutron star cooling, equation 
(21) can be approximated as < Ex >~ 2.7fcT c . 
Using Eqs. (12) and (13) respectively, we have 



-0.1 



< Ex >o.i~ < 



and 



< Ex >o.i~ < 



32.9(f)- 
5.2 x 10 6 (|) 



0.5 



-0.11 



7.49(f)- 
7.33 x 10 4 (|) 



forp_i5>io 2 p 

forp_ 15 <io 2 p 
(34) 

forp_i 5 >i.6P 
forp_i5<i.6P 



( 35 ) 

where P = 10 15 -P-i5 is the period derivative of 
a pulsar in units of s s _1 . In the second case, the 
thermal X-ray come from the bombardment of the 
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relativistic particles from the outer gap. Putting 
Eqs. (28) and (20) into Eq. (29), we have 



f{r, a) 
with 

G{r, a ) -- 



5.2B u 4/7 P w / 21 Rl 0/7 G(r,a) (36) 



cos 9x-y 

W 



r 



a(a) cos a 



Obviously / is a function of r as well as the incli- 
nation angle a in the two cases. 

It is believed that an outer gap start at the null 
charge surface (SI B = 0), which defines the inner 
boundary of the outer gap and the radial distance 
is ri n . From Eq. (32) or Eq. (36), the fractional 
size reaches a minimum at the radius (r in ) of the 
inner boundary, and then increases with radius for 
a given pulsar. Therefore, the fractional size of the 
outer gap at the radius r,„ determine whether or 
not the outer gap exists. If f(r in ,a) > 1, it means 
that the pulsar does not exist any outer gap. In 
other words, a pulsar with f(r in ,a) < 1 should 
have its outer gap and will emit high-energy pho- 
tons produced in the outer gap. As the radius 
increases, the fractional size of the outer gap in- 
creases. For a pulsar with f(r c ,a) < 1, the outer 
gap will extend from r^„ to r c . However, the pul- 
sar with f(r c ,a) > 1 will stop at some radius n,, 
in which /(/■(,, a) = 1. In order to explain the 
average properties of high-energy photon emission 
from the outer gap, we assume that high-energy 
emission at a average radius < r > represents 
the typical emission of high-energy photons from 
a pulsar. The average radius is given by 



< r >= 



f^f(r,a)rdr 

•Cr* f( r ' a ) dr 



(38) 



where r max — min(r c ,rb). The average gap size 
is approximated as /(< r >,P,B) by substituting 
< r > into Eq. (32) and Eq. (36) and in general 
a function of P, B and a because r c is a function 
of P and a, and rb are function of P, B and a. 

Generally, the inclination angles for the pulsars 
are not known well, so we consider the changes of 
the fractional size with the inclination angle. Let 
us define the following function: 

f(<r>,P,B)= V (a,P,B)f (P,B) (39) 



where f a (P, B) = 5.5P 26 / 21 S 12 4/7 is the fractional 
size of outer gap by ignoring the effect of inclina- 
tion angle (Zhang & Cheng 1997). So the effect 
of inclination angle should exhibit in the function 
r](a, P, B). In Fig. 1, we consider the variation of 
outer gap size with magnetic inclination angle for 
various P and B. In panel A of Fig. 1, the upper 
(lower) solid line is the self-consistent outer gap 
f model (the cooling X-ray outer gap model) with 
(37)-P = 0.1s and the upper (lower) dashed line is the 
self-consistent outer gap model (the cooling X-ray 
outer gap model) with P = 0.3s respectively. The 
magnetic field is chosen to be 3 x 10 12 G. In panel 
B of Fig. 1, the upper (lower) solid line is the 
self-consistent outer gap model (the cooling X-ray 
outer gap model) with B\2 = 3.0 and the upper 
(lower) dashed line is the self-consistent outer gap 
model (the cooling X-ray outer gap model) with 
B\i = 1.5 respectively. The period is chosen to be 
0.2s. We can see that r\ varies about a factor of 
2 for the self-consistent model but becomes rather 
constant for the cooling X-ray outer gap model. 

3. High-Energy Gamma-ray Luminosity 

EGRET has observed six pulsars with high 
confidence and three possible radio pulsars to 
emit high-energy gamma-rays above 100 MeV 
(see Thompson 2001 for a review). One of these 
possible radio pulsars is a millisecond pulsar 
PSR J0218+4232 (Kuipcr et al. 2000). There- 
fore, there could be eight young pulsars to emit 
high-energy gamma-rays observed by EGRET. 
It should be noted that PSR B1509-58 is also a 
gamma-ray pulsar. But it is seen only up to 10 
MeV by COMPTEL (Kuiper et al. 1999) and not 
above 100 MeV by EGRET. The observed gamma- 
ray luminosity of a pulsar is L° bs = 4:ird 2 (F 1 , 
where F 1 is the observed gamma-ray flux, d is the 
distance to the pulsar and ( is the gamma-ray 
beaming fraction (0 < £ < 1), which is defined as 
the ratio of the beaming solid angle to An. Two 
parameters are highly uncertain: the distance d 
and the beaming fraction (. It is commonly as- 
sumed that C = 1/47T in estimating the observed 
gamma-ray luminosity. However, the 7-ray beam- 
ing fraction should be different for various 7-ray 
pulsars, which is a function of the magnetic incli- 
nation angle as well as the size of the outer gap. 
Some approximation forms of beaming fraction 
have been given (Yadigaroglu & Romani 1995; Ro- 



7 



mani 1996; Zhang, Zhang & Cheng 2000). How 
accurate of these approximation forms are not 
known. Furthermore, we want to emphasize that 
in addition the beaming fraction the distance es- 
timate also affects the observed gamma-ray lumi- 
nosity strongly. For example, the recent measure- 
ment shows that the distance to the Vela pulsar 
is 294^5Q pc (Cavareo et al. 2001), which is less 
than previous value (500 pc) derived from radio 
observation. Hence the gamma-ray luminosity of 
the Vela could be a factor of 2 lower than that 
given in Thompson ct al.(2001). In this paper, 
we just want to see how the inclination affects the 
gamma-ray luminosity. For simplification, we use 
a common assumption of ( = 1/Att in order to 
compare with the observed data given by Thomp- 
son et al.(2001). 

Because high-energy gamma-rays are mainly 
produced from the outer gap in our model, we 
will compare our expected gamma-ray luminosi- 
ties with those of the gamma-ray pulsars which 
emit gamma-rays observed by EGRET. In our 
model, gamma-ray luminosity for each gamma-ray 
pulsar depends on period, magnetic field and the 
magnetic inclination. However, the magnetic in- 
clination angles are not known well. Once the av- 
erage fractional size of the outer gap for a pulsar is 
estimated, the gamma-ray luminosity can be ap- 
proximated as 



L-y « / 3 (< r >)L sd 



(40) 



where L s d = Air 2 IP jP z is the spin-down luminos 
ity of the pulsar and I — 10 45 g cm 2 . Using Eq 
(4) and putting Eq. (36) into Eq. (40), we have 



with 



L 



7,0 



L jfiV 3 {a,P,B) 



1.36 x 10 33 B% 7 P~ 2/7 



(41) 



(42) 



It is obviously that Eq. (42) is the same as that 
given by Zhang & Cheng (1997). 

In principle, we can estimate 7-ray luminosity 
for a pulsar with a known inclination angle. How- 
ever, it is difficult to estimate 7-ray luminosity 
for all canonical pulsars because only the inclina- 
tion angles of a few of pulsars are known. There- 
fore, we find out statistically the relation between 
7-ray luminosity and pulsar spin-down power for 
the canonical pulsars using Monte Carlo method. 



The details of this Monte Carlo method is given 
by Cheng & Zhang (1998) and Zhang, Zhang & 
Cheng (2000). We use the following assumptions 
for generating the Galactic pulsar population: 

(i) The pulsars are born at a rate (Nns ~ (1 — 
2) per century) with spin periods of Pq = 10 
ms. 

(ii) The initial position for each pulsar is es- 
timated from the distributions p z {z) = 
(l/z oxp )exp(-|z|/z oxp ) a,ndp R (R) = (a R /R 2 xp ) 
R exp(— R/R oxp ), where z is the distance 
from the Galactic plane, R is the distance 
from the Galactic center, z exp — 75 pc, 
a R = [1 - e -»W*«p(l + fl max /#ex P )]-\ 

Rexp = 4.5 kpc and i? max = 20 kpc (Paczyn- 
ski 1990; Sturncr & Dormer 1996). 

(iii) The initial magnetic fields are distributed as 
a Gaussian in logB with a mean log B = 
12.52 and a dispersion gb = 0.35. We ignore 
any field decay for these rotation-powered 
pulsars. 

(iv) The initial velocity of each pulsar is the 
vector sum of the circular rotation velocity 
at the birth location and a random veloc- 
ity from the supernova explosion (Paczyn- 
ski 1990; Cheng & Zhang 1998), the circu- 
lar velocity is determined by Galactic grav- 
itational potential and the random velocity 
are distributed as a Maxwellian distribution 
with dispersion of three dimensional veloc- 
ity cry = \/3 x 100 km/s (Lorimcr, Bailes & 
Harrison 1997). 

(v) A random distribution of magnetic inclina- 
tion angles is used (e.g. Biggs 1990). How- 
ever, the values of a subject to two con- 
straints. First the photon energy given in 
Eq. (28) cannot higher than the total po- 
tential drop of the outer gap. Second n n 
in Eq. (9) must be larger than the stellar 
radius. 



The pulsar period at time t can be estimated 



by 



P(t) 



16ir 2 R% s B 



3/c 3 



1/2 



(43) 
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where, Rns is the neutron star radius and I is 
the neutron star moment of inertia. The period 
derivative (P) can be determined by 

PP = (8tt 2 R% s / '31 c 3 )B 2 . (44) 

Furthermore, the pulsar position at time t is de- 
termined following its motion in the Galactic grav- 
itational potential. Using the equations given by 
Paczynski (1990) for given initial velocity, the or- 
bit integrations are performed by using the 4th or- 
der Runge Kutta method with variable time step 
(Press et al. 1992) on the variables R, Vr, z, Vz 
and <j). Then the sky position and the distance of 
the simulated pulsar can be calculated. 

We now consider the observational selection ef- 
fects. First, we consider radio selection effects in 
order to generate a pulsar population detectable 
at the radio band: the pulsar must satisfy that 
its radio flux is greater than the radio survey flux 
threshold and its broadened pulse width is less 
than the rotation period (e.g. Sturner & Dermer 
1996). The pulsar which satisfies 

Lioo/d 2 > S min (45) 

is considered to be a radio-detectable pulsar, 
where L400 is the radio luminosity at 400 MHz 
and d is the distance to the pulsar. The radio 
beaming fraction can be expressed as (Emmering 
& Chevalier 1989) 

f r (w) = (1 — cosu) + (tt/2 — oj) sinw , (46) 

where u = 6°. 2 x P' 1 / 2 (e.g. Biggs, 1990) is the 
half-angle of the radio emission cone. Then, fol- 
lowing Emmering & Chevalier (1989), a sample 
pulsar with a given period P is chosen in one out 
of fr(P)^ 1 cases using the Monte Carlo method. 
Second, we consider the 7-ray selection effects. 
According to Zhang, Zhang and Cheng (2000), we 
use 

S y {> lOOMcV) > 1.2 x 10~ 10 erg crn^s- 1 (47) 

as the minimum detectable 7-ray energy flux, 
which corresponds roughly to the faintest sources 
with {TS) 1 ' 2 > 5. 

Using above method, we can generate a 7-ray 
pulsar population. In Fig. 2, we show the rela- 
tion between 7-ray luminosity and the spin-down 
power for the simulated 7-ray pulsar population by 



using this new outer gap model, where the birth 
rate of the neutron stars is assumed to be 1/200 yr 
and the pulsar ages are limited to be not greater 
than 10 7 yr. In this figure, we show two pulsar 
populations, one is the population which radio se- 
lection effects of the pulsar with outer gaps are 
taken into account (i.e. radio-loud 7-ray pulsars 
see shaded circles in Fig. 2), and the other popu- 
lation is the radio-quiet 7-ray pulsars (see open 
circles in Fig. 2) . 

There are several interesting features in this fig- 
ure. First there is a rather sharp boundary on the 
left of the population. In fact this boundary is 
given by 

L 7 = L sd . (48) 

This relation/boundary results from the fact that 
7-ray pulsars terminate at /(< r >) = 1 and 
Tmax = H in Eq. (38). Similarly the second rough 
boundary appearing at the bottom of the popu- 
lation is caused by /(< r >) = 1 and r max = 
r c = tl- Zhang & Cheng (1997) has estimated 
the fractional gap size by assuming that the typi- 
cal distance from the gap to the star is <~ rj,, they 
obtained / = 5.5B u 4/r P 26 / 21 . When we substi- 
tute this relation into L sc i, we obtain 

L c ; d lt = 1.5 x lO^P^^rgs- 1 . (49) 

It is important to note that L 7 = L s d as L s d < 

T crit 
L sd ■ 

In such a pulsar population, the best fit gives 
following relation 

log L 7 w 20.42 + 0.38 log L sd (50) 

for the pulsars with outer gaps. When taking the 
radio selection effects, the slope between L 7 and 
L s d becomes flatter, which is log£ 7 w 23.25 + 
0.30 log L s d- We will consider the statistics of 
the pulsars with outer gaps in the rest of this 
section. In order to show the important effect 
of the magnetic inclination angle on 7-ray lu- 
minosity, we show the change of log 10 (L 7 /i 7i o) 
with log 10 L 7 in Fig. 3. The best fit indicates 
that log 10 (L 7 /£ %0 ) oc 0.91og 10 L 7 . For compari- 
son, we show the result given by Zhang & Cheng 
(1997), which is independent with the inclination 
angle. 

It is very interesting that the slope between L 1 
and L s d will become steeper when the minimum 
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detectable 7-ray energy flux decreases. For ex- 
ample, we obtain a 7-ray pulsar population using 
GLAST threshold (1.8 x 1(T 12 erg cm- 2 s" 1 ), the 
best fit gives 

log £ 7 w 17.45 + 0.46 log L sd . (51) 

1 II 

That is L 7 cx LJ d . In this pulsar population, we 
have log 10 (-L 7 /L 7 ,o) oc 0.92 log L 1 . 

For comparison, we show our results given by 
Eqs. (50) and (51) and the observed high-energy 
photon luminosities of eight 7-ray pulsars in Fig. 
4. In this figure, the observed data are taken from 
Thompson (2001), which derived from detected 
fluxes above 1 eV assuming the solid angle of pho- 
ton beaming of 1 sr. From Fig. 5, it can be seen 
that both our results given by Eqs. (50) and (51) 
are consistent with the observed data. For the re- 
sult fitting the simulated 7-ray pulsar population 
with EGRET threshold, the expected slope (0.38) 
is flatter than the observed one (0.46). 

Finally we would like to point out that this new 
outer gap model predicts many more weak 7-ray 
pulsars whose 7-ray luminosities can be as low as 
10 32 erg s^ 1 (cf. Fig. 2) than previous outer gap 
models (Zhang & Cheng 1997; Cheng & Zhang 
1998). The main reason for existing these weak 7- 
ray pulsars is that instead of using the outer gap 
size at half of light cylinder radius to determine 
if the outer gap can exist, the new model allows 
the outer gap to survive when the fractional size 
of outer gap at the null charge surface for a given 
pulsar is less than unity. This prediction allows 
the age of 7-ray pulsars to extend to a few million 
years old. In fact, more detail Monte Carlo simu- 
lation results show that most of these weak 7-ray 
pulsars have ages between 0.3-3 million years old 
and many of them are located at higher galactic 
latitude. Most of these weak 7-ray pulsars their 
radio beam could miss the Earth and they could 
contribute to the unidentified 7-ray sources in high 
galactic latitude (Cheng et al. 2003). These pre- 
dictions can be verified by GLAST. 

4. The Death lines of Pulsars with Outer 
Gaps 

We now consider the condition which the outer 
gap of a pulsar exists. For an outer gap, the inner 
boundary is estimated as the null charge surface 
where the magnetic field lines are perpendicular to 



the rotation axis. From Eq. (36) or (32), f(r,a) 
reaches minimum at r = r{ n . In other words, if 
the fractional size of the outer gap at r in is larger 
than unity, then the outer gap would not exist. 
Therefore, we can estimate the death lines of the 
pulsars with outer gaps by using f{r ini a) = 1. 
It should be noted that G pc (ri n , a) in Eq. (32) 
or G(ri n ,a) in Eq. (36) is only the function of a 
because r in /B,L only depends on a (see Eq. (9)). 
For the case of X-rays produced by the neutron 
star cooling and polar cap heating, we obtain from 
f(r in ,a) = 1 

j o p_ ( 3.1 logP — Ai(a) foriogP<i3+io g p 
— \ If- log P - A 2 (a) forio g p>i3+iogp 

(52) 

for the temperature given by Zhang & Harding 

(2000) , where A x (a) = 12.87- 3.16 logG pc (a) and 
A 2 (a) = 12.92- ^\ogG pc (a); and 

^ p _ ( 3.1 log P — Bi (a) forio g p<i4.8+iogp 
°^ \ |logP — B 2 (a) forio g p>i4.8+iogp 

(53) 

for the temperature used by Hibschman & Arons 

(2001) , where Bi{a) = 13.46- 3.16 logG pc (a) and 
B 2 (a) = 13.95 - 21ogG pc (a). For the case of 

X-rays produced by the bombardment of the re- 
turning particles from the outer gap, we have 

10 7 
logP = —logP- 13.02+ -log G(a) . (54) 

Obviously, the death lines depend on the magnetic 
inclination angle. 

Although the magnetic inclination angle of each 
pulsar has not been determined well, the distri- 
bution of the magnetic inclination angle can be 
estimated from the statistical polarization study 
of the radio pulsars. It was generally believed 
that the parent distribution of the magnetic in- 
clinations satisfy an uniform distribution (Gunn 
& Ostriker 1970; Gil & Han 1996). However, re- 
cent study by using polarization data of the radio 
pulsars indicate that the parent distribution of the 
magnetic inclinations satisfies a cosine-like distri- 
bution (Tauris & Manchester 1998). Therefore, 
we estimate the average value of /(rj„, a) in these 
two possible parent distributions of the magnetic 
inclinations, i.e. 

< G(a) >= [ G(a)U(a)da/ [ U{a)da , (55) 
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logP 
and 



{ |1< 



where U(a) represents the distribution of the in- 
clination angles. We consider uniform and cosine- 
like distributions of the inclination angles respec- 
tively. For the uniform distribution, we have 
< G pc {a) >« 0.32 and G(a) = 0.38. Therefore 
Eqs. (52), (53) and (54) become 

. p_( 3.1 logP - 14.43 forio g p<i3+iogP 
g \ f log P- 13.77 foriogP>i3+io g p , 

(56) 

logP— 15.02 forio g p<i4.8+iogp 
logP— 14.94 forio g p>i4.8+iogp 

(57) 

logP = ^ log P- 14.60 . (58) 

For cosine distribution, < G pc (a) >m 0.43 and 
G(a) w 0.49, Eqs. (52), (53) and (54) become 

p_ f 3.1 logP - 14.03 forio g p<i3+iogp 
g \ f log P- 13.55 forio g p>i3+iogp , 

(59) 

^ p _ ( 3.1 logP - 14.62 forio g p<i4.8+iogp 
X |logP— 14.68 forio g p>i4.80+iogp 

(60) 



and 



logP = ^ log P- 14.20 . (61) 



Based on the original outer gap model (CHR 
I; CHR II), Chen & Ruderman (1993) have con- 
sidered the death lines of Crab-like and Vela-like 
pulsars respectively. For the Crab-like pulsars, the 
death line is (see Eq. (25) of their paper) 



logP = 41ogP-7 . 



(62) 



For Vela-like pulsars, they introduced a param- 
eter £ to describe the expected variation of the 
magnetic field within the outer gap and assumed 
£ = 1/2. In this case, the death line is (see Eq. 
(27) of their paper) 



log P = 3.8 log P- 10.2 . 



(63) 



According to Zhang & Cheng (1997), we have from 

/, « 2.83 X l -4p20/21p-2/7 = ! 

logP = ^ log P- 12.42 . (64) 

In Fig. 5, we show the death lines of the pul- 
sars with outer gaps for two possible distributions 



of the magnetic inclination angles. In this case, we 
assume that X-rays are produced by the neutron 
star cooling and polar cap heating. Because the 
temperature used by Hibschman & Arons (2001) 
is higher than that used by Zhang & Harding 
(2000), the average energy of X-rays for the for- 
mer is greater than that for the latter. Therefore, 
the death lines obtained by using the tempera- 
ture of Hibschman & Arons (2001) are lower than 
those by using the temperature of Zhang & Hard- 
ing (2000). 

In Fig. 6, we show the death lines of the pulsars 
with outer gaps, in which X-rays are produced by 
the returning relativistic particles from the outer 
gap, and the outer gap is self-sustained. For com- 
parison, we also plot the death lines given by Chen 
& Ruderman (1993), and the death line derived 
from the model of Zhang & Cheng (1997). It can 
be seen that our model predicts that much more 
radio pulsars have self-sustained outer gaps com- 
pared to those given by Chen & Ruderman (1993) 
as well as by Zhang & Cheng (1997). In our es- 
timate of the death lines of the pulsars with self- 
sustained outer gaps, we require that /(r in ) = 1. 
This condition is reasonable. It is believed that 
an outer gap can develop along the null charge 
surface (Cheng, Ruderman & Sutherland 1976) or 
along the last closed field line (CHR I; CHR II). 
Therefore, a self-sustained outer gap exists if the 
fractional size of the outer gap at r in is not greater 
than unity. 

5. Conclusion and Discussion 

After taking the geometry of the dipole mag- 
netic field, we have given a revised version of the 
outer gap model given by Zhang & Cheng (1997). 
In the revised outer gap model, the fractional size 
of the outer gap is not only the function of period 
and magnetic field of the neutron star, but also 
the function of the radial distance (r) to the neu- 
tron star and the magnetic inclination angle (a). 
In other words, the fractional size of the outer gap 
has a form of f(r,a) = /»(P, Bi 2 )Gi(r,a), where 
/o(P, P12) is only the function of pulsar period 
and surface magnetic field, Gi{r,a) changes with 
radial distance to the neutron star and the mag- 
netic inclination angle and subscript i represents 
the X-ray field considered. The fractional size of 
the outer gap is given by Eq. (32) in the X-ray 



11 



field which is produced by the neutron star cool- 
ing and polar cap heating and by Eq. (36) in the 
X-ray field which is produced by outer gap heat- 
ing. In this model, the fractional size of the outer 
gap has a minimum at the inner boundary for a 
given pulsar, increases with the radial distance 
along the last open field lines and then reaches 
its outer boundary where /(r^a) = 1. In other 
words, the outer gaps of some pulsars do not start 
from the null charge surface to the light cylinder. 
We have shown the changes of the fractional size 
of the outer gap with the 

inclination angle (see Fig.l). Further, we have 
found that the outer gaps of relative young pulsars 
such Vela can extend from the null charge sur- 
face to the light cylinder for any inclination angle 
, however, the outer gaps of some pulsars like as 
Geminga cannot extend to the light cylinder for a 
larger inclination angle (say 75°). 

In order to describe the average properties of 
high-energy radiation from a 7-ray pulsar, we have 
defined an average radial distance < r > (see Eq. 
(38)). Using Monte Carlo method described by 
Cheng & Zhang (1998) (also see Zhang, Zhang 
& Cheng 2000), we simulated two populations 
of 7-ray pulsars whose energy fluxes are greater 
than EGRET threshold and GLAST threshold. 
In these simulations, we only consider the case of 
pulsars with self-sustained outer gap and uniform 
distribution of the inclination angles. We have 
plotted the change of L 7 with L s d (see Fig. 2). 
We also indicated the variation of L 7 /L 7 .o with 
L 7 in Fig. 3, which show the importance of the 
inclination angle on L 1 . In the model of Zhang 
& Cheng (1997), L 7 is independent on the incli- 
nation angle. Fitting the simulated results, we 
found that L 1 oc L°J S for EGRET threshold and 
i 7 oc L°J 6 for GLAST threshold. Compared with 
the observed data given by Thompson (2001), our 
simulated results are reasonable (see Fig. 4). In 
fact, current distance of the Vela pulsar (Cavareo, 
ct al. 2001) is less than that used by Thompson 
(2001) about a factor of two, which reduce the lu- 
minosity about a factor of four. 

In Fig. 4, we note that Geminga is not within 
the simulated population. Zhang & Cheng (2001) 
have used a three-dimension outer gap model to 
explain the phase-resolved spectra of Geminga. 
They discovered that in order to fit the observed 
7-ray data the solid angle Aft must be near 5sr. 



In Fig. 4, Ail = Isr is used for all observed 
7-ray pulsars. If the larger solid angle is used, 
it brings Geminga within the simulated popula- 
tion. It should be pointed out that Yadigaroglu 
& Romani (1995) have studied the effect of 7-ray 
beaming in their outer gap model (also see Romani 
1996). Zhang et al (2000) also gave a approximate 
expression of the 7-ray beaming fraction, which 
will be applied to estimate 7-ray fluxes in our new 
model. 

It is interesting to point out that both polar gap 
models and outer gap models predict L 1 = L s d for 
low spin-down power pulsars. But the position of 
the break occurs at 5 x 10 33 P- 1/2 er#/s for the 
polar gap model (Harding et al. 2002) whereas 
the outer gap models predict a higher position at 
1.5 x 10 34 P 1 / 3 er(/s _1 . For higher spin-down power 
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pulsars, the polar gap models predict L 7 = L sd , 
whereas the outer gap cannot give precise predic- 
tion because the model L 7 also depends on a less 
well-known parameter a. The statistical predic- 
tions of the outer gap model give L 7 = L s sd , where 
6 depends on the properties of 7-ray detector. For 
example, S = 0.38 for the EGRET and <5 = 0.46 
for GLAST. 

According to our model, the fractional size of 
outer gap at the null charge surface for a given 
pulsar (f(r in ,a)) reaches a minimum. Therefore, 
the outer gap should exist only if f(ri n ,a) < 1. 
Averaging f(ri n ,a) on two possible (uniform and 
cosine) distributions of the magnetic inclination 
angles respectively, we have obtained the death 
lines of the pulsars with outer gaps in the two pos- 
sible X-ray fields and the comparison them with 
the observed data in Figs. 5 and 6. 

Compared with the death line derived from the 
outer gap model of Zhang & Cheng (1997), the 
revised model predict that more pulsars will have 
their outer gaps and then emit high-energy pho- 
tons. 

We would like to make the conclusion as follows. 
Our results indicate that (i) the intrinsic param- 
eters for explaining the observed 7-ray properties 
of rotation-powered pulsars are the magnetic incli- 
nation angle, period and magnetic field. We have 
obtained a very concrete functional form of the 
prediction of gamma-ray luminosity which only 
depends on these intrinsic parameters, the incli- 
nation angle could be known if the radio data is 
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sufficiently good; (ii) the conversion efficiency for 
7 known gamma-ray pulsars, that is rather scat- 
tered , can be explained in our revised model using 
the Monte Carlo method. Although our estima- 
tion of the outer gap size cannot precisely predict 
the thickness of individual pulsar when the mag- 
netic inclination angle is poorly known, it is a rea- 
sonable estimation of the statistical properties of 
gamma-ray pulsars using the statistical method; 

(iii) Unlike those 7 observed gamma-ray pulsars, 
mature pulsars (ages 0.3-3 million years) can also 
be gamma-ray pulsars and their efficiency is insen- 
sitive to the inclination angle. Most importantly, 
their gamma-ray luminosity is proportional to the 
spin-down power, which can be tested by GLAST; 

(iv) The mean cut-off age of gamma-ray pulsars is 
increased by a factor of 3, therefore older gamma- 
ray pulsars (up to about 3 million years old) can 
move up to higher galactic latitude. Some uniden- 
tified EGRET gamma-ray sources could be mature 
pulsars with ages between 0.3-3 million years old 
pulsars predicted by this model. In fact, Parkes 
Observatory survey has discovered a large num- 
ber of radio pulsars on the error boxes of EGRET 
unidentified gamma-ray point sources (Torres et 
al. 2003); and (v) The previous works on death 
line based on the outer gap model (e.g. Chen 
and Rudcrman 1993) did not give a detail physical 
reason rather they used a phenomcnological ap- 
proach. Here we proposed a new criterion based on 
concrete physical reason for the death line, which 
predicts more gamma-ray pulsars. (vi)Our model 
predicts many more weak 7-ray pulsars with age 
between 0.3 3 million years old than previous 
outer gap models and they satisfy L 7 oc L s d, which 
can be verified by GLAST. 

Finally, we would like to remark that the 
gamma-ray luminosity formulae developed in this 
paper may not be able to explain the gamma-ray 
luminosity of individual pulsar. Two important 
factors, i.e. distance uncertainty and beaming 
fraction, which play crucial roles in determining 
the gamma-ray luminosity, have not been con- 
sidered here. For example, the luminosity ra- 
tio between PSR B1055-52 and Geminga, which 
have same spin-down power of 3 x 10 34 er<? s _1 , 
is about a factor of 6 (Thompson et al. 2001). 
But the model prediction (cf. Fig. 2) is no more 
than a factor of 3. In fact, the estimate of the 
distance of PSR B1055-52 is very uncertain, its 



distance can change from ~ 1.5 kpc (for exam- 
ple, Thompson et al. 2001) to a small value of 
^500 pc (see, for example, Ogelman & Finley 
1993; Combi et al. 1997; Torres, Butt, Camilo 
2001; Hirotani & Shibata 2002), which makes 
the ratio change from <~ 10 to <~ 3. Recent dis- 
tance estimate for PSR B1055-52 from the dis- 
persion measure is ~ 0.72 kpc (see http: //rsd- 
www.nrl.navy.mil/7213/lazio/ne_model, also see 
Mignani, DcLuca & Caraveo 2003). However, it 
is well known that 25% error is common in de- 
termining the dispersion measure (MaLaughlin & 
Cordes 2000), which gives a factor of 2 uncertain- 
ties in the luminosity. Also the distance estimate 
of Geminga is known to have at least 25% uncer- 
tainties (Caraveo et al. 1996), which give another 
factor of 2 uncertainties in gamma-ray luminos- 
ity estimate. Other uncertainty results from the 
gamma-ray beaming fraction. In principle, both 
distance estimate and beaming fraction will be 
obtained more accurate. Then our model still 
predicts that the difference in the magnetic in- 
clination angle can still cause a large difference 
in gamma-ray luminosity for pulsars with same 
spin-down power. 
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Fig. 1. — Variation of the fractional size of the 
outer gap with the magnetic inclination angle 
for some typical pulsar parameters. Panel (A): 
7](P, B, a) versus a for a given magnetic field of 
3 x 10 12 G. The upper (lower) solid line is the 
self-consistent outer gap model (the cooling X-ray 
outer gap model) with P = OAs and the upper 
(lower) dashed line is the self-consistent outer gap 
model (the cooling X-ray outer gap model) with 
P = 0.3s respectively. Panel (B): r/(P,B,a) ver- 
sus a for a given period of 0.2 s. The upper 
(lower) solid line is the self-consistent outer gap 
model (the cooling X-ray outer gap model) with 
B12 = 3.0 and the upper (lower) dashed line is the 
self-consistent outer gap model (the cooling X-ray 
outer gap model) with Bi 2 = 1.5 respectively. 



Fig. 2. — The change of 7-ray luminosity (_L 7 with 
the spin-down power (L st j) in the 7-ray pulsar pop- 
ulation predicted by our outer gap model. In our 
simulation, we have used the EGRET threshold as 
the minimum detectable 7-ray energy flux. Open 
circles and shaped circles are the model radio-quiet 
and radio-loud 7-ray pulsars respectively and the 
solid line is the best fit for all 7-ray pulsars with 
outer gaps. Shaded line is the best fit for the radio- 
loud 7-ray pulsars with outer gaps. 
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Fig. 3. — Lj/Ljfi versus L 1 in the 7-ray pulsar 
population predicted by our outer gap model. In 
our simulation, we have used the EGRET thresh- 
old as the minimum detectable 7-ray energy flux. 
Open circles are the expected data and solid line 
is the best fit. For comparison, we show the result 
given by Zhang & Cheng (1997) as a dashed line. 
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Fig. 5. — Death lines of the pulsars with outer 
gaps. It is assumed that X-rays are produced 
by the neutron star cooling and polar cap heat- 
ing. Two cases of the inclination angle distribu- 
tions are considered: (i)an uniform distribution 
and (ii) a cosine distribution. Dotted and solid 
lines are presented the death lines given by Eqs. 
(56) and (59). Dashed and dot-dashed lines repre- 
sents the death lines given by Eqs. (57) and (60). 
The observed data are taken from see website 
http://www.atnf.csiro.au/research/catalogue/) 



Fig. 4. — 7-ray luminosity versus the spin-down 
power. Solid circles with error bar are the ob- 
served data given by Thompson et al. (2001), solid 
and dashed lines are our results given by Eqs. (50) 
and (51) respectively. 
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Fig. 6. — Death lines of the pulsars 

with self-sustained outer gaps. The ob- 
served data are taken from see website 
http://www.atnf.csiro.au/research/catalogue/. 
Lines (1) and (2) are given by Eqs. (62) and 
(63) respectively (Chen & Ruderman 1993). Line 
(3) is the death line (Eq. (64)) predicted by 
Zhang & Cheng (1997). Lines (4) and (5) are 
the death lines of our model for the uniform and 
the cosine distributions of the inclination angles, 
respectively. 
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